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We study a quantum model with non-isotropic two-dimensional oscillator 
potential but with additional quadratic interaction x\x<i with imaginary cou- 
pling constant. It is shown, that for a specific connection between coupling 
constant and oscillator frequences, the model is not amenable to a conven- 
tional separation of variables. The property of shape invariance allows to find 
analytically all eigenfunctions and the spectrum is found to be equidistant. 
It is shown that the Hamiltonian is non-diagonalizable, and the resolution 
of the identity must include also the corresponding associated functions. 
These functions are constructed explicitly, and their properties are investi- 
gated. The problem of R— separation of variables in two-dimensional systems 
is discussed. 

1. Introduction. 

During last years, starting from a pioneering paper of C. Bender and S.Boettcher p], there 

is a growing interest to investigate Quantum Mechanics with non-Hermitian Hamiltonians 

(see also [2]) consistently. It was shown that under definite assumptions the spectrum of 

such Hamiltonians is real and a modified scalar product which provides unitary evolution 

a E-mail: cannata@bo.infn.it 
b E-mail: m.ioffe@pobox.spbu.ru 
c E-mail: cutaisi@yahoo.com 



M. V. Ioffe 2 R D. N. Nishnianidze 23 @ 



1 



can be built for some models. For comprehensive reviews, see papers [3], [I] and references 
therein. 

With few exceptions [5], [6], [7], the analysis concerned one-dimensional Quantum Me- 
chanics. In particular, most results were obtained for a wide class of models, with unbroken 
PT— invariance PQ, [8], [9], [10] . It can be considered as a modern generalization of conven- 
tional Quantum Mechanics to a non-Hermitian one. 

In turn, the notion of the pseudo-Hermiticity: 



with r] a Hermitian invertible operator, allowed to define a more general class of non- 
Hermitian systems with physically acceptable properties of energy spectra. The most sys- 
tematic investigation of pseudo-Hermiticity has been performed by A.Mostafazadeh [11] (see 
also [12], [H]). A suitable description of Hilbert space for such systems is given in terms 
of biorthogonal basis, which consists of the eigenstates |\l/ n ) and \if? n ) of H and iP, corre- 
spondingly. 

It was found that some systems with complex potentials are naturally described by Hamil- 
tonians which are not diagonalizable. They correspond to the systems whose biorthogonal 
basis elements do not provide complete basis in Hilbert space. In such a case, one has 
to add the so-called associated functions to complete the basis, and Hamiltonian becomes 
block-diagonal with some number of Jordan blocks of standard structure on its diagonal. 

In Section 2 we formulate the two-dimensional model with complex potential having the 
form of second order polynomial in xi,X2- Usually, such model is solved easily by means of 
linear transformation of coordinates with subsequent separation of variables, maybe complex. 
This procedure was described, for example, in [6] . But two peculiar cases of such polynomial 
potentials with special relation between constants are beyond this scheme: they are not 
amenable to separation of variables. Just such model is studied in Section 2. It is solved 
exactly: the whole energy spectrum and corresponding wave functions are found analytically. 
Instead of by separation of variables, which is impossible here, the model is solved by means 
of shape invariance, a powerful method introduced in the framework of SUSY Quantum 
Mechanics. In Section 3 we investigate the properties of the constructed wave functions. We 
show that they do not realize a resolution of identity, i.e. they do not form complete basis, i.e. 
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the Hamiltonian is not diagonalizable [14]. The corresponding associated functions are also 
built analytically in this Section, and their properties are studied in detail. In Section 4, we 
discuss the conventional procedure of separation of variables in two-dimensional Schrodinger 
equation both with real and with complex potentials. For the first case, the old results of 
Eisenhart [15] are reproduced, and for the second case, we prove that the model of previous 
Sections does not allow the most general (nonlinear) algorithm of R— separation of variables 

m- 

2. Two-dimensional complex oscillator. 

Let us consider the two-dimensional model with complex oscillator Hamiltonian: 

H = -A (2) + V(x) = -d\ -dl + u\x\ + tujx 2 2 + 2igx 1 x 2 . (2) 
Performing the linear complex transformation of variables X\,x 2 

2 

Xi = ^ ^ a ijVjy (3) 
J'=l 

where are complex elements of matrix A, one may try to separate variables in the 
Schrodinger equation: 

H^f(x) = E^{x). (4) 

It is necessary to obtain in diagonal form both the Laplacian and the quadratic potential, 
though in complex variables yi. As one can check, this is possible for generic values of 
parameters Wj, (7 in (J2j), with two exclusions. Indeed this is impossible, if the coupling constant 
is: 

2g = ±(u 2 1 -u 2 2 ), (5) 

when the Jacobian of ([3]) vanishes. Just this situation will be considered below in this paper, 
and for definiteness we will choose the minus sign above. 

We will use the complex variables z = x\ + ix 2 , z = X\ — ix 2 = z*, for which 

H = -Ad z d- Z + \ 2 zz + gz 2 - 2 A 2 = cu 2 + u\ > 0; A > 0. (6) 
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One can check that the Hamiltonian fl6]) obeys the following property: 

HA + = A + {H + 2X); HA~ = A~ {H - 2 A) , (7) 
with operators A^ defined as: 

A ± = d z ^ \z- (8) 

The equation (jTj) realizes the particular case of shape invariance p2] , the property appeared 
on the first time in the framework of one-dimensional SUSY Quantum Mechanics [18] with 
real potentials. Shape invariance of the Hamiltonian H(x; a), which depends on a parameter 
a means that this Hamiltonian satisfies the intertwining relations with some operators Q^: 

H(x;a)Q + = Q + H{x;a) + R(a); a = 5(a); R(a) > 0; (9) 
Q-H(x;a) = H(x; a)Q~ + R(a). (10) 

This property provides a very elegant method to solve the Schrodinger equation algebraically: 
practically, all known one- dimensional exactly- solvable models are shape-invariant. Recently 
this property was generalized to the two-dimensional Quantum Mechanics [19], where it gives 
usually a quasi-exact-solvability of the model (i.e. analytical construction of a part of wave 
functions). For the case of two-dimensional Morse potential with specific values of parame- 
ters, shape invariance helped to find the whole spectrum and corresponding eigenfunctions 

The Eq.(J7|) corresponds to the simplest variant of shape invariance: a(a) = a and R(a) = 
2A. This case was investigated in one- dimensional Quantum Mechanics |2T] . and the same 
idea must work also in two-dimensional case. Since the intertwined Hamiltonians in ([9]) 
coincide, such Hamiltonians were called self-isospectral, and for one- dimensional case self- 
isospectrality leads to an equidistant (oscillator-like) character of the spectrum. In our 
present case of two-dimensional complex shape invariant potential, the property (EJ) provides 
also the oscillator-like spectrum of H : 

E n = 2X(n + l). (11) 

We deal here with unusual Quantum Mechanics - with complex potential and, even more, 
with non-diagonalizable Hamiltonian (see the next Section). Nevertheless, the absence of 
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singularities in fl2]) ensures that nothing like well known "fall to the center" phenomena 
[22] is possible here (the formal proof can be found in Appendix). The normalizable bound 
state wave functions will be exponentially decreasing at infinity, having no singularities. The 
corresponding spectrum is bounded from below, the ground state with energy Eq is denoted 
as ^o t o(x) (it will be clear below, why we use two indices for enumeration of ty). The excited 
levels correspond to n — 1, 2, ... in fITT]) . 

From the second intertwining relation ([7]) it follows that \l/o,o must be a zero mode of 
A- : 

A-* ,oOM) = 0. (12) 

Otherwise, the Hamiltonian H would have the lower level (Eq — 2A), contrary to our assump- 
tion on \l/o,o above. The solution of zero mode equation ( JT2l) can be found in general form. 
Indeed, it can be written in terms of a new function of z only: 

*o,o(z, z) = exp (—~zz)4f (z), 

and the Schrodinger equation (j3J) leads to the first order differential equation: 

2Xz%(z) = (E -gz 2 -2\)V Q (z). 

Its solution 

^o(^) = (^) (i?0 - 2A)/2A exp[-(^ 2 )] 

has to be a single- valued function on a plane, i.e. to be 2tt— periodic in polar angle <p leading 
therefore to the spectrum (fTTlh Formally, infinitely many functions solve (112)): 

*n,o(2, z) = c n! oZ n exp(-azz -bz 2 ); n = 0,1,2...,; a = -; b = ^ ( 13 ) 

with normalization constants c nj o- The actual ground bound state corresponds to n = in 
( TT3l) . and it has energy 2 A. 

All excited states can be built by the standard algebraic procedure of shape invariance: 

* n (z, z) = (A + ) n V 0fi (z, z); n = 1, 2, ... (14) 

with exactly the same energies as in (fTTj) . It is easy to check, that expressions in ( !T3|) and 
( IT4|) coincide for all n > 0, i.e. \l/ n (-2, z) = ^^0(2;, f ), being the zero modes of operator A~ . 
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One can also prove that assuming the existence of any eigenstate different from ( Fl4|) we 
would obtain the level with real part lying below E Q , in contradiction with the definition of 
the ground state having the lowest real part of the energy. Thus, the whole spectrum of the 
system (J2j) is exactly known - ( flTT) . The corresponding wave functions fy n ^(z,z) in (fl"3l) are 
known analytically as well, but additional investigation of their properties is required. It 
will be performed in the next Section. 



3. Non-diagonalizability of the Hamiltonian. 
3.1. General scheme. 



The Quantum Mechanics with non-Hermitian Hamiltonians needs a suitable modification 
of the scalar product and resolution of identity to make the model self-consistent [3], [I]. In 
general, if the non-Hermitian Hamiltonian H commutes with some antilinear operator B, 
this may be used for definition of a new scalar product between arbitrary elements of Hilbert 
space as follows: 

= J (Btt)*. (15) 

In the ordinary Quantum Mechanics with real potentials just the customary complex con- 
jugation plays the role of such operator B, and the Hamiltonian H is Hermitian under such 
scalar product (\E , |<&). 

In the case of a general antilinear operator B, the Hamiltonian H obeys Hermiticity but 
with a scalar product (|T5l) . Symbolically: 

((V\H\$)) = J(B^)H<$> = J(HBy)$ = J(BHV)$= ((HV\$)), (16) 

where double integration by parts and vanishing of off-integral terms were used to move H 
to the left under the integral. In the two-dimensional case, the Ostrogradsly-Gauss theorem 
(a two-dimensional analog of two integrations by parts) allows to move H as well, due to 
exponential decreasing of all wave functions on large contour. Many one- dimensional models 
with non-Hermitian Hamiltonians obey the so-called PT— symmetry pQ, [3], citebenderOl, 
[9], where the role of B is served by the antilinear symmetry operator PT of simultaneous 
time and coordinate reflections. In more general situation, non-Hermitian Hamiltonians may 
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be pseudo-Hermitian [3] (the definition in operatorial form was given in (0Q)). Then the opera- 
tor B can be chosen as B = r/T, and the scalar product (|T5|) coincides with r\— scalar-product 
commonly used in the literature. The pseudo-Hermitian Hamiltonians ([I]) are Hermitian be- 
ing considered under this scalar product. 

In our case of two-dimensional coordinate space, such antilinear operator, keeping H in 
02]) invariant, can be chosen either as P{F or as P2T, where Pi means X\ -H- —X\, and P2 
means £2 -H- — x-i. Let us choose the second option for definiteness. The wave functions 
\l/ n >0 (z, z) are simultaneously the eigenfunctions of P 2 T with unique eigenvalue +1. For such 
choice of the operator B, the scalar product ((\I/|<I?)) becomes an integral over the product 
J instead of the J in the ordinary Quantum Mechanics. 

Now the important property of the wave functions under a new scalar product has to be 
investigated - are the corresponding norms positively definite or not. The norms of the basic 
states \l/ nj o can be calculated explicitly: 



Thus, only the ground bound state \&o,o is normalizable state with a positive norm. All ex- 
cited wave functions have zero norms (are "self-orthogonal"), prohibiting the usual resolution 
of identity in terms of complete set of eigenfunctions of H. Such situation was investigated 
in one- dimensional Quantum Mechanics with complex potentials in a series of papers [14J: 
the zero norm of wave function signals that one deals with a non-diagonalizable Hamilto- 
nian. Then it is necessary to build the so-called associated functions which participate in a 
resolution of identity and complete the basis. 

The adequate formalism for diagonalizable non-Hermitian Hamiltonians is the so-called 
biorthogonal basis in the Hilbert space [3], [11]. This basis includes two types of states 
l^n), \^ n ), such that: 




(17) 



H\* n ) = E n \V n }; #t|* n ) = E*\* n ); (* n |* m ) = (* m |*„) = 



nm 



(18) 



7 



with decompositions: 

J = E|* n )(* n |; H = EE n \* n )(* n \. (19) 
In coordinate representation, one can take ^ n (x) = \l/*(x), so that the scalar product is: 

(*„|* m ) = J y n (x)^ m (x)d 2 X = ((*n|*m)) = Sum- (20) 

This formalism has been discussed in [J], [11], [13], for a particular explicit calculations see 
also [23] (in one- dimensional case) and [7] (in two dimensions). 

The formalism becomes more complicated [14J in the case of quantum systems with non- 
diagonalizable non-Hermitian Hamiltonians, like our (j2]). For such systems, the basis (|T8l) is 
not complete, and the decompositions (|T9l) do not work. In order to improve the situation, 
every self-orthogonal wave function \l/ n) o, n > 1 with zero norm must be accompanied with 
a set of p n — 1 associated functions "$? n ,k, k = l,2,...,p n — 1. It must be clear now, why 
notations with two indices of wave functions were introduced above. By definition, these 
functions obey: 

(H - E n )V nik = *„, fc _i; k = l,2,...,p n -l, (21) 

where all functions are supposed to be normalizable, and the operator H, when applied to 
these functions, preserves this property. Here we restrict ourselves for simplicity with the 
case when each self-orthogonal eigenfunction ^^O; n — 1 ; 2,... is accompanied by only one 
set of associated functions ^ n ,k, k — 1,2, ...,p n — 1. 

Similarly to the scheme of previous paragraph, the partner eigenfunctions ^ n ,o also are 
accompanied by their associated functions ^ n> k, k = 1,2, ...,p n — 1. Practically, it is conve- 
nient to numerate the functions identifying them with as follows: 

V n , Pn - k -i = K,k fc = 0,l,2,...p n -l. (22) 

With these notations, according to the general formalism which was illustrated in detail for 
some one- dimensional models [H] , the scalar product (|20|) in the extended biorthogonal basis 
must be: 

«*„,*!*„,,,» = (*^> = / •m^W* = W 

k = 0, 1, ...,p n -I] I = 0, 1, ...,p m - 1. 
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Correspondingly, the generalized decompositions must be: 

oo pn— 1 
n=0 fc=0 

OO p n — 1 OO p n — 2 

= ^^K|^n, fc ))((^,pn- fc -l|+^^l^))((^,pn-fc-2|. (25) 

n=0 fc=0 n=0 k=0 

The Hamiltonian is clearly non-diagonal, but block-diagonal. Each block - Jordan cell of 
standard form (see ff25l) ) - has dimensionality p n . 

3.2. The specific model: Non-Hermitian two-dimensional oscillator. 



The eigenf unctions \l/ n o were found analytically in the previous Section (see Eq. fflB"]) ). All 
these functions for n > 1 were shown to be self-orthogonal, and therefore, some associated 
functions must be properly taken into account. In this Subsection we are going to investigate 
the properties of these functions, and in particular, to check the relations ( 1231) and to find 
the values of p n . 

First of all, we will prove that scalar products (1231) vanish for different energy levels 
E n , E m , i.e. that ({^ n ,k\^ m.i)) = for n ^ m. The proof is by induction in k, I with 
essential use of pseudo-Hermiticity of H. Indeed, as in ordinary Quantum Mechanics: 

= ((^, |#|^ m ,o» - «jy*^ |*m,0» = (E m ~ E n )((y nfi \Vm,o)), 

i.e. wave functions with different energies are orthogonal. Analogously, because of definition 
= ((^ nfi \(H-E m )\^ m ^))-(((H-E m )^ nfi \^ mA )} = 

= ((^n,o|*m,o)) - (E n - E m )({9f n>0 \9 m>1 )), 

and the scalar products between wave functions and first associated functions for different 
E n , E m also vanish: 

«*n,o|*m,l» = «*n,l|* m ,o)> = 0. 

The procedure can be continued further leading to orthogonality of all functions with different 
n, m. 



9 



Now we have to consider the wave functions and associated functions, which belong to 
the same value n. For the first associated function ^ ny \, the defining equation (12T]) can be 
solved explicitly in a general form with two arbitrary constants: 



0"n, \ZZ — a n> i — - — Z + C n: iZ + —{C n fi — )z m z 

exp [— (— zz + ^—z 2 )]. 
2 4A n 



One of integration constants is defined by physical requirement for wave functions to be 
single- valued in the plane: 

Ac„0 C n Q , , 

^ = " ^' (26) 



leading to: 



<M*> z) = [^zz n - 1 - (n ifi 1) ^ z- 2 + c n>1 z n ] exp [-{azz + bz% (27) 
8o lbab 



The second integration constant c n> i reflects the obvious fact (see (}2"Tl)). that ^ n> i is defined 
up to a solution of homogeneous Schrodinger equation c n ^ n $. This additional term must 
be defined by a suitable "normalization", i.e. by conditions ( l23l) . One can check, that the 
normalization conditions ( |23l) fix uniquely not only c n> i, but all higher additional terms as 
well. 

To find the norm of ^ n ,i and some other scalar products below, we shall calculate a class 
of two-dimensional integrals of the form: 



Ijv.m = J z N z M exp [— 2{azz + bz 2 )]dzdz 



(2f 



with positive constants a, b and integer A/", M. These integrals vanish for odd values of (N + 
M) due to antisymmetry under a space reflection (xi,x 2 ) — >■ — (a;i,x 2 ), i.e. ijy,M = for 
(AT + M) = 2s + 1. For even values of (N + M), we start from the basic integral [23]: 

/(a, 6, c) = J exp [-2(az* + fez 2 + cz 2 )]dzdz = thT 1 ; 5 = 2 ^ (a 2 -46c). (29) 



10 



Then, the required In,m can be calculated by suitable differentiations of I(a, b, c) = 1(8) : 
hn,2{n+k) = j \zz) 2n z 2k exp [— 2(azz + bz 2 + cz 2 )]dzdz\, 



■ c=0 



= \{-\da) 2n {-\d b ) k I{5)}\ c=Q = 0; k > (30) 
h{n+k),2n = j \zz) 2n z 2k exp [-2(azz + bz 2 + cz 2 )]dzdz\ c=0 = 

i(-\d a ) 2n (~d c ) k I(5)] lc=0 = n(-l) k 2-( 2n+1 \2k + l)\(2k + l) 2n b k a~ (2k+2n+1 ^; (31) 

hn+i,2(n+k)+i = J {zz) 2n+l z 2k exp [-2(azz + bz 2 + cz 2 )]dzdz\ c=0 = 

= [{-\da) 2n+ \-\d b ) k m]\c=o = 0;k>0 (32) 

h(n+k)+i,2n+i = J {zz) 2n+1 z 2k exp [-2(azz + bz 2 + cz 2 )}dzdz\ c=0 = 

= \(-\da) 2n+ \-\d c ) k I(5)}\ c=0 = n(-l) k 2-( 2n+2 \2k + l)\(2k + l) 2n+l b k a^ 2k+2n+ ^3) 

In particular, it is clear that In,m = for M > N. The listed integrals allow to check that 
for the first excited level n = 1 functions ^1,0, ^1,1, in addition to self-orthogonality of ^1,0? 
satisfy (up to normalization factors) the relations (1231) with the value p± = 2 : 

^1 1^1 dzdz = 1; / ^1 i^ndzdz = 0, 



where one has to take c^i = in (1271) and = 32a 2 b/ir in ([TBI . For the general value 

n > 2, the norm of the first associated functions j (^> n ^) 2 dzdz = Const ■ 5 n ^- Comparing 
with the orthogonality relations ( [231 . this signals that p 2 = 3, and vice versa, p n 7^ 3 for 
n > 2. 

Calculation of scalar products for higher order associated functions ^ n ,k with k > 2 will 
be performed by an alternative method, where the main role is played by the last associated 
functions ^ niPn -i. From the definition ( J2TT) . the following expression for ^ n .k in terms of 
^n )Pn -i can be easily derived: 

^ n , fc = (H — E^-^r^- k = 0, 1, ...,p n - 1. (34) 

The dimension p n — 1 of Jordan cell and the highest associated function ^ n ^ n _i must be 
derived by solving the equation: 

(H - E n y^ niPn ^ = 0, (35) 
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by subsequent calculation of \l/ n ,fc along and finally, by checking the required scalar 

products (1231) . 

We shall prove now that p n = n + 1 and the following solution of Eq. ()35l) 



ty n! n = c n Q n exp[—(azz + bz 2 )}, Q n = (az + bz) n , n>0, (36) 

satisfy all necessary conditions above (constants of normalizations c n will be defined below). 

In order to prove these statements, one has to check straightforwardly the following 
relations: 

(H - E m ) ■ ^0,0 = ^o,o • 4 \-d z d- z + a(zd z + zd- z - m) + 2bzd z ) = ^o,o ■ An, (37) 



(D m Q n ) = 4ay— bn(n — l)f2 n _2 + (n — m)Vt n + 2bnzQ n -\ J , (38) 

An ■ & = z j D m - m 3 ~ X Qz + *ajz>. (39) 
By the method of mathematical induction, it can be derived from (1371) -(1391). that: 

D k n VL n = (2ab) k (n - (2k - i - 1)\ otf h%. 

i=0 V / 2k-i 



'n-(2fc-i)j 



where: 

(a) k = a(a + !)...(« + k - 1); a?> = {-2)\ a« = 2 2 *, a£> <fc = (~ 2 )^ ~ ^ 
and the coefficients a) satisfy the system of equations: 



2(i - k)af ] = {i + l)a[% < i < k] 
2(2aS-af)) = a f +1 ), 1 < i < k. 



Thus, 

(H-E n ) k ^ n ^ n = c n (2ab) k exp[-(azz + bz 2 )}Y,^ k) (n-(2k-t-l) ) ^.(afc-i), (40) 

i=0 ^ 



2k-i 
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and for k = n, only one coefficient of z l does not vanish: % = k = n. Therefore, 



(H - E n r* n:n = c n (8ab) n n\z n exp [-(azz + bz 2 )} = Cn ^ a ° } U ^ nfi . (41) 

c n,0 

The action of one more operator (H — E n ) onto this equation leads just to required Eg. ([55]) . 
The norm of ^ n ,n is calculated using the explicit expression fl3T)|) and the integrals flBTJ]) - 

((^ n ,n\^n,n)) = J dzdz(az + bz) 2n exp [-2(azz + bz 2 + cz 2 )] c=Q = 

= J dzdz(-^) n (a 2 d c + 2abd a + b 2 d b ) n exp [-2{azz + bz 2 + cz 2 )} c=0 = 

= (-\Vi(a 2 d c + 2abd a + b 2 d b ) n I(6)] c=0 = 0. 

The action of the operator (ad a + bd b ) onto the vanishing integral above leads to the useful 
expressions for arbitrary integer n, k : 

dzdz(az + bz) n z k exp [-2(azz + bz 2 )] = 5 nk ——. (42) 

Equations ( 140]) for k < n, together with Eqs.( l34p and ( 142)) . allow to analyze other scalar 
products of associated functions with the same energy Indeed, since according to ( 140]) \l/ nin _fc 
is a linear combination are combinations of constants entering (140p ): 

k 

^n,n-k = PiZ l VLri-2k+i^ Ofli 

i=0 

the required scalar products are: 

((^n,n-k,^n,n-k')) = / dzdz ^ PiPj?* 3 ^n~2k+i^n-2k' +j exp [-2 (dZZ + fe 2 )] = 

^ i=0 j=0 

= f dzdzj2J2^ i+i ( az+b ^ 2n ~ 2(k+k ' )+ ^ 

J i=0 j=Q 

Therefore, due to Eq.( )42|) . the following choice of c nj0 in ffl3|) and c n in f )36|) : 

c nj0 2 2a(16a6) n 

° n = (8a6)«n!' C "'° = * 
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provides normalized scalar products 



((*n,fc,*n,fc'» =8 k ,n-k', (43) 

as it should be for p n = n + 1. Summarizing this Subsection, we proved that the dimension of 
Jordan cell corresponding to energy E n depends on n, namely, p n = n+1, and all constructed 
associated functions, after suitable normalization, provide the necessary scalar products f l23l) . 

4. Non-separability of variables. 



It was already noted in Section 2, that the Schrodinger equation with (complex) potential 
of second order in xi, X2, just in the case ([5]), does not allow an ordinary separation of 
variables by means of linear transformation of coordinates. A more general question will be 
studied in this Section: whether any nonlinear transformation can provide the separation. 
Actually, we are interested in opportunity to perform the so-called R— separation of variables 
[16], which is explored, for example, in three-dimensional problems with central forces. In 
two-dimensional context, R— separation means that coordinates Xi, x 2 can be mapped (and 
the mapping is invertible) to new (complex, in general) variables: 

qx = Qi(xi, x 2 ), q 2 = Q 2 (x x ,x 2 )\ xi = X 1 (q 1 ,q 2 ), x 2 = X 2 (q 1 ,q 2 ), (44) 



so that the Hamiltonian takes the form: 
1 



H 



(45) 



n(gi) + r 2 (q 2 ) 

with arbitrary functions r», /ij, In such a case, the problem splits onto two one- dimensional 
problems. The separation of variables in two- and three-dimensional Schrodinger equation 
with real potentials was investigated by L.P.Eisenhart |15] , where the exhaustive lists of 
corresponding coordinate systems qi(x) and potentials V(x) were found: eleven systems 
in three dimensions and three systems in two-dimensional case. Because of complexity 
of potential, we are interested here in generalization of these results on two-dimensional 
R— separation: both potentials V(x) and the new coordinates qi,q 2 may be complex. 



14 



The change of variables (jHJ) in the kinetic part of H = —(d^ + d X2 ) + V(xi,x 2 ) gives 
the following conditions: 

{d Xl Qi){d xl Q 2 ) + {d X2 Q x ){d X2 Q 2 ) = (46) 
(d Xl Qi) 2 + (d, 2 Qi) 2 = (d Xl Q 2 ) 2 + (d x . 2 Q 2 ) 2 = ] (47) 

r l \Ql) + T 2{C[2) 

which lead, in particular, to relation: 

d xl Q 2 = d X2 Q x . (48) 

Actually, the opposite sign in the r.h.s. of ( 1481) is also possible, but the final results will be 
the same. The general solution of (T4"8"|) is expressed in terms of an arbitrary complex function 
G: 

Qx = d X2 G(x 1 ,x 2 ); Q 2 = d Xl G(xi,x 2 ). (49) 
After substitution of (SSj) back into fT46l) . 

(d Xl d X2 G(xi, x 2 )) ■ (d 2 Xl + d 2 X2 )G(x u x 2 ) = 0, (50) 

we have two options: 

d xl d X2 G(x 1 ,x 2 ) = 0; (51) 

(d 2 Xl +d 2 X2 )G(x u x 2 ) = 0. (52) 

For the first option,the variable q± depends on X\ only, and analogously, q 2 on x 2 . Then 
from (1471) . 

2 / \ 2 



^iQi(^i) = d X2 Q 2 (x 2 ) 



1 



|a:i=JCi(ffi) V 2 /|iea=X2(«a) T l(qi) + T 2 (q 2 ) 

for which only non-interesting solutions exist: Ti(qi) = Const, and correspondingly, Qi(xi 
are linear complex functions of X{. 

The second option (|52|) gives more interesting general solution: 

G(xi, x 2 ) = m(z) + n(z); z = x\ + ix 2 ; z = x x — ix 2 

q\ = Qi(xi,x 2 ) = m\z) + n'(z); q 2 = Q 2 (xi,x 2 ) = i(m'(z) - n'{z)); 

z = z (Q-)\ z = Z(q+) = (z{q~)j ; q+ = qi + iq 2 ] q- = q x - %q 2 
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in terms of new functions m,n, Z, Z . In order to investigate, whether variables (71,(72 allow 
separation of variables in the Schrodinger equation, we substitute this solution in the kinetic 
part of Hamiltonian: 

(d 2 Xl + d 2 X2 ) = Am(z)n(z)(d 2 qi + d 2 ^)\ z=z{q _)-- z =z{ q +y 
The separation of variables ( I45p in Laplacian is possible only if 

(4m(z)n(z) |2=z(g _ );5=%+) ) _1 = 71(91) + r 2 (g 2 ), 

i.e., if: 

d qi d q2 {m-\Z(q_))n-\Z(q + )) =0. 
Thus, one obtains two ordinary differential equations: 

K^f= ( ^f = 7 2 ; 1 = const. (53) 
n {q+) m {q-) 

After straightforward calculations, the case 7 = leads to the following solutions: 

z — ot\q 2 _ + (3iq~ + «i; z = a 2 q\ + (3 2 q + + k 2) (54) 

where cti, fy, Ki are constants. Here z = x\ + ix 2 , z = x± — ix 2 are still mutually conjugate, 
although q + , q_ are, in general, not necessarily conjugate. If we are interested, similarly to 
[j~5] . only in real variables qi, and therefore q + = g* , then the constants in (l54"j) are mutually 
conjugate: ot\ = a 2 = a, /?i = /3 2 = «i = ^2 = ft- In this case for a 7^ by means of 
a suitable constant shifts of qi and Xj, without loss of generality, one can made (3 = k = 0. 
Choosing also the special scale, namely, a = 1/2, we obtain: 

2 2 

g 2 2 = ±xi + (x 2 + x 2 ) 1/2 ; xi = gl 2 ; x 2 = q x q 2 , (55) 

and 

r{qi) = q-] + r 2 (g 2 ) = 2(x 2 + x 2 ) 1/2 . (56) 

The potential takes the form: 

V( Xl , x 2 ) = {x\ + x 2 )- 1 / 2 [/(x 1 + (x 2 + x 2 ) 1 / 2 ) + g(- Xl + {x\ + x 2 ) 1 / 2 )]. (57) 
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Expressions (155]) - (1571) coincide exactly with the case III of Eisenhart [15] . 

For «i = «2 = in ([Ml), the transformation x — > q describes linear transformations 
to complex coordinates qi,q2, which were mentioned in Section 2. For complex qi,q2 the 
general form ( l54j) is allowed. For example, a\ = 0, a 2 7^ among others. 

The case 7 7^ gives from ([53]) : 



z = o x exp (79-) - 5i exp (—79-); -2 = cr 2 exp (7?-) - 5 2 exp (-79-). (58) 

For the case of real q\,q 2 , the constant 7 is real, and <j\ = a 2 = er, #i = #2 = 8. Two options 
for ( [58]) must be considered separately. If 5 = 0, by means of rotations in (x±, X2) — plane one 
obtains: 

.Ti = cr exp (7^! ) cos 7g 2 ; x 2 = -a exp (7^) sin ^q 2 \ 
T(9i) + r 2 (g 2 ) = tf 2 exp (2791) 

T/(x 1 ,x 2 )= /(X2/3:i) 2 + ^ + ^ ) , (59) 



coinciding with the case I of Eisenhart 

The second option with real a = —5 = a/2 gives exactly the case II of Eisenhart: 

Xi — a cosh (751) cos (792); x 2 = ~ a sinh (791) sin (7^2); 
7^ct^ 

7i(gi) + r 2 (q 2 ) = —— [cosh (2791) - cos (27g 2 )]; 

After the analysis above, it is clear that R— separation of variables for polynomial poten- 
tial in 02]), if possible at all, would belong to the option ( 151]) . But explicit substitution of 
([54]) into shows that this expression is not reducible to the form (H5[) . i.e. the system ([2]) 
is not amenable to separation of variables. 

5. Conclusions. 



It is interesting to compare the non-Hermitian model ([2]), ([5]), which does not allow 
for separation of variables, with the same model, but without condition ([5]), i.e. with the 



17 



model of [6j. It is clear that our model corresponds to merging of pairs of mutually complex 
conjugate energy levels in [6], since the restriction fl5]) just leads to vanishing imaginary parts 
of energy eigenvalues. Thus, the condition fl5]) on coupling constant violates diagonalizability 
of the Hamiltonian Eq.(3) in [6] and also turns complex energy levels to the real axis. The 
special remark concerns the dimensionality of Jordan cells: p n — (n + 1) coincides with the 
degeneracy of the corresponding levels E n in the model ([2]), for g = and U\ = U2, i.e. in 
the model of isotropic real two-dimensional oscillator. It would be interesting to investigate 
further the classical integrals of motion and most importantly quantum symmetry operators 
for the diagonalizable and non-diagonalizable cases, both being solvable, and the interplay 
with separability of variables. 
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Al. Symmetry operator. 

It follows from the shape invariance property, i.e. from intertwining relations (I7|), that 
the operator R = A + A~ commutes with Hamiltonian H. This fact allows to choose the wave 
functions of H such that they are simultaneously the eigenf unctions of R. Let us solve the 
corresponding pair of equations: 



where eigenvalues of R are denoted as — A 2 r 2 /4 with arbitrary (complex) constant r. The 
equation f lSTj) has two independent solutions: 




(61) 
(62) 



{-Ad z d- Z + X 2 zz + gz 2 )^ = 



^ (1) (^, z) = c {1 \z) exp(zf(z)), ¥ 2 \z, z) = c {2 \z) exp(-zf(z)), 



(63) 
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where c^-' 2 >{z) are arbitrary functions, and f(z) = \\J z 2 — r 2 /2. Substitution of (163]) into 
( 162]) leads to first order homogeneous equations for c^' 2 \z) : 

-4(c (1) /y + (gz 2 -E)c {l \z) = 0, (64) 

4(c (2) /)' + (gz 2 -E)cW(z) = 0. (65) 

They are solvable explicitly, in a general form: 

c^iz) = ^exp( f 9 *~ E dz), c&>(z) = ^exp(- / ^ ~ E dz). (66) 



/ 4/ " ^ / ^ J 4/ 

It is easy to check that ^( 2 \z,z) has an exponentially decreasing asymptotics at infinity. 
The apparent singularity of c^(z) at zeros of f(z) can be compensated for r = just for 
the positive values of E, thus justifying the choice of positive n in ( JTTj) . For these values of 
energy, the functions ^( 2 \z,z) coincide with wave functions (1T31) in Section 2. 
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